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Abstract: I study the zero-temperature phase transition between superffuid and insulating 
ground states of the Bose-Hubbard model in a random chemical potential and at large integer 
average number of particles per site. Duality transformation maps the pure Bose-Hubbard 
model onto the sine-Gordon theory in one dimension (ID), and onto the three dimensional 
Higgs electrodynamics in two dimensions (2D). In ID the random chemical potential in 
the dual theory couples to the space derivative of the dual field, and appears as a random 
magnetic field along the imaginary time direction in 2D. I show that the transition from the 
superfluid state both in ID and 2D is always controlled by the random critical point. This 
arises due to a coupling constant in the dual theory with replicas which becomes generated 
at large distances by the random chemical potential, and represents a relevant perturbation 
at the pure superfluid-Mott insulator fixed point. At large distances the dual theory in 
ID becomes equivalent to the Haldane's macroscopic representation of disordered quantum 
fluid, where the generated term is identified with the random backscattering. In 2D the 
generated coupling corresponds to the random mass of the complex field which represents 
vortex loops. I calculate the critical exponents at the superfluid-Bose glass fixed point in 2D 
to be v = 1.38 and z = 1.93, and the universal conductivity at the transition o c = 0.25el/h, 
using the one-loop field-theoretic renormalization group in fixed dimension. 



1 Introduction 



In recent years, a new paradigm in physics of condensed matter has emerged under the name 
of superconductor-insulator transition. Among the quantum many-particle systems which 
are believed to exhibit this type of transition between their ground states at zero temperature 
are 4 ife in random media |l|, thin superconducting films |2| and Josephson junction arrays 
||. At T = 0, as some parameter in Hamiltonian is varied, such a system is expected to show 
either zero or infinite dc linear resistance, with the transition caused by purely quantum 
fluctuations which qualitatively alter the many-body ground state. Experimentally, the 
superconductor-insulator transition is manifested as a change from continuously increasing 
to sharply decreasing resistance of the system as temperature is lowered, in accordance with 
the notion that at T = the ground state is either an insulator or a superconductor. 

It has been argued ||], that as a model for superconductor-insulator transitions 
it suffices to consider a Hamiltonian for interacting bosons in random external potential. 
The basic assumption behind the idea is that the transition corresponds to the onset of 
phase coherence of the already preformed Cooper pairs, and not to the formation of the 
pairs themselves. The ground state of the bosons can then be either superfluid (SF) or 
insulating; the insulating state may arise due to repulsive interactions, in which case it is 
an incompressible Mott insulator (MI) with a gap, or from a combination of interaction and 
localization effects, which may result in the formation of the gapless Bose glass phase (BG). 
Particularly interesting is the SI transition in two dimensions (2D). At finite temperature 2D 
are the lower critical dimension for superconductivity with the superconducting transition 
being of Kosterlitz-Thouless type ||. Also, for non- interacting particles 2D are the lower 
critical dimension for Anderson localization, with all single particle states weakly localized jrj. 
Arguably, the main source of interest in 2D superconductor-insulator transitions comes from 
the suggestion || that right at the quantum (T = 0) critical point dc conductivity should 
be finite and universal, determined only by the universality class of the transition. This 
exotic possibility of normal but universal diffusion results from a delicate balance between 
localization and superfluid fluctuations, and is unique for 2D. 

Although the experimental support for universality of the conductivity at the transition 
is still rather weak []2|, the calculation of this and other critical quantities for different 
universality classes of dirty-boson systems in 2D is a fundamental and unsolved problem. 
Without disorder, the transition in the system of interacting bosons on a lattice is always 
between Mott insulator and superfluid phases. With integer average number of particles per 
site, MI-SF transition is in the universality class of the classical 3D XY model ||, ||. The 
critical exponents in this case are well known to be z — 1, v « 0.667, and the universal 
conductivity has been calculated by variety of methods, including 1/N and the Monte Carlo 



calculations ||10|| , and the e-expansion ||11|| . With disorder present, the critical behavior at 
the superconductor-insulator transition in 2D is much less understood, due to apparent non- 
existence of the upper critical dimension where the theory would have a critical point at a 
weak coupling |J, [12]. In fact, even the very nature of the transition at commensurate boson 
densities has been a matter of debate: simple arguments || suggest that at weak disorder 
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the transition into the superfluid phase proceeds from the Bose glass phase, while most of the 
numerical studies see a direct MI-SF transition Jl3| , (T^j . In absence of a controlled analytical 



approach to the problem with disorder, most of our knowledge of the critical behavior at 
the BG-SF transition originates from the numerical studies |15|], [16|], p!7fl . On the analytical 



side, large-N expansion [JIB] , real space methods |19| and strong coupling expansion [^(J have 
all been brought to bear. While these calculations give valuable informations on the phase 
diagram of the system, the critical exponents are difficult to extract and less reliable. In 
particular, the universal conductivity at the BG-SF transition has hitherto been calculated 
only numerically. 

In the present paper I study the critical behavior at T = superconductor-insulator 
transition in ID and 2D systems of interacting lattice bosons in Gaussian random chemical 
potential and at commensurate densities within dual description of the problem |21| which 
focuses on topological defects in the many-body ground state. This approach, in which 
one looks at the destruction of the long-range order by proliferation of the defects, has 
proven invaluable for several problems in the past, the Kosterlitz-Thouless transition being 
a prime example. Let me first describe the main results for the simpler ID system. By 
duality the disordered Bose-Hubbard model maps onto the sine-Gordon theory in which the 
random chemical potential couples linearly to the space derivative of the dual field. The 
elimination of the fast modes of the dual field performed perturbatively in the strength of 
the lattice potential always leads to exponential suppression of the small periodic potential 
by disorder, so that at some smaller cutoff the lattice potential in the theory scales to 
zero. Most importantly, while renormalizing to zero the periodic potential together with the 
random chemical potential generates an additional disorder term, so that at long length scales 
the dual theory is without the periodic potential term, but contains two different disorder 
terms. This long length scale effective theory is equivalent to the Haldane's macroscopic 
representation [p2]| , || of the disordered quantum fluid in ID, where the generated term arises 
from the backscattering of the random potential. Once it helped generate the important 
backscattering term in the action, the random chemical potential which corresponds to 
forward scattering in the long length scale theory becomes redundant and can be removed 
from the problem by an appropriate change of variables. Thus the ID disordered Bose- 
Hubbard model at large length-scales is described by the Haldane's macroscopic theory for 
dirty bosons with no trace of the lattice potential. The transition from the superfluid state 
is even at infinitesimal disorder controlled by the random fixed point in ID ||23|| , ||, which is 
interpreted as a sign that the transition is between the superfluid and the Bose glass phases. 
I end the analysis of the ID problem with the discussion of the shape of the BG-SF phase 
boundary in the interaction- disorder plane. 

The lesson from the ID example is that the transition from the superfluid state in the 
disordered system is always into the Bose glass, and that it is controlled by the fixed point 
which can be located easily in the space of coupling constants of the dual theory. This 
suggests the same line of attack at the 2D version of the problem. A short account of the 
results in 2D has already been published pi| . The dual theory for 2D dirty bosons at T = 
and with short-range repulsion is considerably more complicated than its ID counterpart, 
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and has the form of 2+1D classical, anisotropic Higgs electrodynamics (Ginzburg-Landau 
superconductor) in random magnetic field, directed and correlated along one dimension. The 
complex Higgs field describes closed vortex- loops in space and imaginary time, which interact 
via long-range forces mediated by the gauge field, whose Maxwell term in the action describes 
the sound excitations in the parent superfluid state. Condensation of the vortex field implies 
destruction of the superfluid by appearance of infinitely large vortex-loops, which is also 
reflected in the induced Higgs mass of the gauge field. At the pure critical point which 
controls the MI-SF transition, weak anisotropy which arises from the underlying quantum 
nature of the problem is argued to be irrelevant. More interestingly, weak disorder in the 
form of random correlated magnetic field turns out to be precisely marginal, and I assume it 
is marginally irrelevant. Similarly as in the ID case, the random chemical potential generates 
additional disorder which can be understood as the random "mass" of the vortex field, which 
is strongly relevant at the pure MI-SF fixed point. The flow at the critical surface eventually 
ends in the attractive random critical point, which I again interpret as the BG-SF fixed 
point. 

The critical exponents and the universal resistivity at the BG-SF critical point in 2D 
are determined to the lowest order in fixed-point values of the coupling constants in the 
dual theory, using the field-theoretic renormalization group in fixed dimension [f2q| . The 
calculations are performed directly in 2+1 dimensions since the duality between the Bose- 
Hubbard Hamiltonian and the Higgs electrodynamics holds only for the specific dimension, 
and I also want to determine the universal conductivity at the critical point, which is again 
finite exclusively in 2+1D. The idea of this method is to express the universal quantities like 
the critical exponents or the universal conductivity as perturbative series in renormalized, 
instead of bare, coupling constants, so that they assume finite values at the transition, 
determined by the critical point of the theory. The lowest order results z = 1.93 and v = 1.38 
obtained this way compare reasonably well with the numerical calculations [|17]] , even though 
I do not have a truly small parameter. The conductivity is obtained using the Kubo formula 
for the response of the disorder field to an external potential. By duality of charges and 
vortices this way one directly calculates the boson resistivity, instead of the conductivity, 
as a series in critical coupling constants. Inverting the final result I obtain the universal 
conductivity a c = 0.25(e^//i) at the quantum SF-BG critical point, in rough agreement 
with the numerical results [[RJ, [17]], and apparently smaller then the value suggested by 
experiments. Subtleties involved in interpretation of the experimental measurements of o c 
are briefly discussed. 

The paper is organized in the following manner. In the next section I define the dirty- 
boson lattice model and use duality transformations to arrive at the lattice theory in density 
representation in general dimension. In section 3 I study the dual theory in ID and discuss 
the renormalization scenario and the phase diagram. In section 4 the dual representation 
in 2D is derived. The one-loop renormalization of the continuum version of the dual theory 
is performed and the critical exponents and the conductivity at the BG-SF transition are 
calculated. In the last section the summary and the discussion of the main conclusions is 
given. Technical details are left for appendices. 
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2 Hamiltonian and duality transformation 



I study the Hamiltonian for a system of charge-e* bosons with on site repulsion U > 0, 
written in number-phase representation (Josephson junction array Hamiltonian) [Bj, ||17| , 
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where index i labels the sites of a ID (2D) lattice, v = x [y = x, y) is a unit vector, 
hi represents a deviation from a large integer average number of bosons per site, hi is a 
Gaussian random chemical potential with hihj = WhSi.j and hi — 0, and 0j is the phase 
variable canonically conjugate to the number of bosons, [4>i,hj] = idij- The Hamiltonian 
(1) is in the same universality class as the Bose-Hubbard model [p6|| , |T7| . I calculate the 
partition function of the system 

Z = Trexp{-(3H} 1 (2) 

P = 1/T, in the basis of states that diagonalize the particle numbers: 1 = J2{ ni } \{ n i}) ({ n i}\i 
[rii] = {ni, ...u Nl }, — oo < rii < oo and integer, where N L is the total number of lattice 
sites. Following the standard procedure for construction of the path-integral representation 
of the partition function [27], one arrives at the expression: 



,7T N L N e JJ 

z = yi j n n ^,« ex p-{^-n n ?a _e H^ n ^~ ( 3 ) 

{rtj iQ } !r i=la=l i,a i,a 

te cos(4>i +0ta - <pi )a ) +i^^, a (n i]a+1 - n i)Q ,)}, 



with the boundary conditions only on the number of bosons n i = n it ^, e = /3/N, and the 
limit e —>■ is assumed. At this point there are two ways to proceed: utilizing the Poisson 
summation formula one may perform the sums over numbers of bosons and arrive at the 
phase representation, from which it follows that without disorder the MI-SF transition at 
T = is in the universality class of the XY model @, ||. Disorder in this formulation enters 
via an imaginary boundary term, reducing significantly the utility of the phase representation 
in that case. Nevertheless, the phase representation leads naturally to the field theory for 
the superfluid order parameter, studied extensively in the literature [12||. Alternatively, one 



can attempt to integrate the phases out to be left with the partition function in number- 
representation. To that end I use the Villain representation of the cosine-term in the last 
equation: 

2 

771 ■ 

exp{et ^ cos(4>i +i> , a - ijCt )} ~ ex P{~ 2et + i .? m w i> ( < ^+ £ '> a _ 4>i,a)}, (4) 

where rn,j ja are one-component (two component) integer valued currents in ID (2D) at the 
sites of the two (three) dimensional lattice labelled by indices (i,a). The integration over 
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the phases now can be performed exactly, and its effect is to produce the constraint on the 
currents and numbers of bosons: 



■ m i>a + V T n^ a = 0, (5) 

where Vf and V r are the lattice gradients (finite differences) in space and imaginary time 
directions, respectively. The last equation may be interpreted as a continuity equation, and 
the two (three) component vector Mj jQ , = (m^ a , rii )a ) represents the conserved bosonic current 
in ID (2D). The Eq. (5) ensures that only closed current loops in space and imaginary time 
contribute to the statistical sum. The partition function now assumes the form 

Z= £ ^V-^exp-^^n^-e^^ + i-^mJJ. (6) 

{M ia } i,a i,a 6 i,a 

Note that the action in the exponent is completely real. As it stands, the last form of the 
partition function is still difficult to study analytically due to discreetness of the variables, 
and is more amenable to numerical methods |17) . It is possible however to further transform 



the partition function into a more telling expression P£j| , [p9 |, [30], as in this context was 
pointed out by Fisher and Lee The form of the transformation differs in ID and 2D, so 



I will address the two cases separately. 

3 One dimension 
3.1 Dual theory 

In ID the constraint V • M = may be resolved by introducing a unique integer scalar 
variable N, so that M = (V T N, —V X N). The theory then takes the form 

Z = lim / rj^, Q exp-{^^(V,A Q ) 2 + ^tE( V -A«) 2 (7) 

-e V hi(V x A i:Ct ) V cos(2vrA i>Q )}, 

y 

i,a a i,a 

where the limit y — > serves to force the real- valued fields {Ai a } to take only integer values. 
Hereafter I will take T = 0. Assuming that softening this constraint will not change the 
universality class of the transition |2"5|], p0||, in the continuum limit a — > 0, e — > the dual 



theory becomes: 

S = K [ dxdr[(d T e(x,r)) 2 + c 2 (d x 6(x,T)) 2 } - I dxdrh(x)d x 6(x J T) 



-v J dxdr cos28(x,t) 



5 



where ttA(x,t) = 6(x,t), c 2 = Uta 2 , K = l/(2tir 2 a), and v may be understood as a (finite) 
strength of the lattice potential with the lattice constant a. Averaging over the Gaussian 
random potential by introducing replicas in standard way [31], one finally arrives at the form 



S = K Y. dxdT[(d T 9 a (x, t)) 2 + c 2 (d x 9 a (x, r)) 2 ] - (9) 

a=l ^ 

N r N . 

w h / dxdrdr'd x 9 a (x,T)d x 9 l3 (x,T f ) — / dxdr cos29 a (x 1 r), 

afi=\ a=l 

where Greek indices now enumerate replicas, and the limit N — > at the end of the cal- 
culation is assumed. If Wh = 0, dual to the Bose-Hubbard model (1) in ID is the familiar 
sine-Gordon theory which has the transition in the XY universality class as K or v are var- 
ied @. At large K the configurations with 9 close to constant dominate the statistical sum, 
which implies a sharp particle number at each site and therefore disordered phases. Notice 
that by simple power counting the disorder coupling Wh is relevant at the pure fixed point 
(located at K*c* — 1/4, v* — 0, see below), so one may expect that with disorder the nature 
of the transition in the system should be changed. I will demonstrate that indeed that is 
the case, although, interestingly, the coupling Wh will ultimately be irrelevant for the critical 
behavior. It is however crucial for bringing about the correct form of the effective theory at 
the large distances, to which I turn next. 



3.2 Renormalization 

I now integrate out the components of the field 9 in Eq. (9) with the momenta A/s < \k\ < A, 
A ~ 1/a and s ~ 1, and with any frequency u, perturbatively in strength of the periodic 
potential v. The calculation is standard, with the only non-trivial step being a matrix 
inversion needed to obtain the propagator for the dual field 9, which in presence of disorder 
is non-trivial in replica indices. The interested reader is refereed to the Appendix A for 
details, and here I only mention that this can be done exactly in the limit — > 0. To the 
lowest order in v the recursion relations are: 

iv P-^-^)« + o(A do) 



d\n(s) v 2Kc 2K 2 c 4 
dK dc 



0{v% (11) 



d\n(s) d\n(s) 

w h + 0(v 2 ), (12) 



dw h ^, 2 



<iln(s) 

and I introduced the dimensionless combinations w^/A — > Wh, v/A 2 — > v and nK — > K. 
At Wh = one may chose the units of lengths and energy so that c = 1, so the Eqs. (10)- 
(12) reduce to the celebrated Kosterlitz recursion relations |§ with K playing the role of 
temperature for the equivalent classical 2D XY model. The crucial new feature brought by 
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disorder is the presence of the u^-term in the recursion relation (10) with a negative sign. 
Physically this is exactly what one expects: at large distances the random chemical potential 
washes away the periodic lattice potential. In fact, since Wh is a relevant perturbation at 
the pure fixed point (at v* = 0, K* = 1/4), we see that it will make a weak enough periodic 
potential exponentially irrelevant at large distances for any K. If this would be the end of 
story we would come to a somewhat paradoxial conclusion that disorder always turns Mott 
insulator into a superfluid. This would follow from the observation that once v = in (9) 
one can eliminate the uvterm altogether by a shift in field 8(x, r) ||23|| , no matter how large 
Wh might have become at that scale. This would be incorrect however, for the following 
reason: to the second order in v a new disorder term in the action (9) becomes generated 
during the mode elimination 

JV 



~D( S ) J2 I dxdrdr' cos 2 (0 a (x, r) -6 p (x,r')), (13) 

a, /3=1 J 

with 

v 2 w h ln(s) 

D{S) = 4KW ' (14) 
where v and Wh are now left dimensionful. The reader will recognize the term (13) as 
precisely the one appearing in Haldane's course-grained representation for the problem |[22|| , 
0, and which arises there from the components of the random potential with the wave 
vector k ~ 2irpo and may be understood as the backscattering term. The w^-term would 
correspond to purely forward scattering in the large length scale theory in ID. Once the 
coupling constant D has become generated and small v has scaled to zero we may freely 
shift the dual variables to completely eliminate the w/j-term from the problem. At a smaller 
cutoff A' the effective theory therefore acquires the familiar form: 

S A , = K(A') J2 I dxdr[(d T 6(x, r)f + c(A') 2 (dJ(x, r)) 2 ] (15) 
-D{A')J2 I dxdrdr' cos 2(6 a {x, r) -6 p {x,t')), 

with the effect of integration over the fast modes absorbed into the effective values of the 
coupling constants. For weak initial v the value of -D(A') at the scale at which v vanishes 
will be small, so one may monitor the change of the effective couplings in (15) during further 
mode elimination perturbatively in D. The recursion relations for this case are well known 
|2§: 

dD 1 , . . 

dHsj = (3 -K )D > (16) 

= ^ (17) 

where again I put c = 1. The transition from the superfluid state is thus always ultimately 
controlled by the random fixed point at D* = and K* = 1/3, at which the correlation 
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length exponent v = oo and the dynamical exponent is z — 1. This is in accord with the ref. 
H where the same conclusion was reached using only slightly different reasoning. Similar 
arguments have also recently been advanced in ||32|| . The random fixed point has a natural 
interpretation as the BG-SF critical point, although the nature of the insulating phase can 
not be directly inferred from the recursion relations at small D. Although the exponents 
have the same values as at the MI-SF fixed point in ID, it is a different critical point: for 
instance, the superfluid two-point correlation function at the transition decays algebraically 
at large distances with the power 1/3, instead of 1/4 at the MI-SF transition [23|. 



3.3 The phase diagram 

Although the precise phase diagram for the ID Bose-Hubbard model will depend on the 
microscopic details, one can still understand some of its general features on the basis of the 
above renormalization scenario. Let me first assume that the periodic potential is weak, that 
c = 1, and think of the phases in K — Wh plane. In these units K = ^JU/t/2TT 2 . At zero 
disorder, for K < K c « 1/4 the system is in the superfluid phase, and for K > K c it is a 
Mott insulator. Now take the point with K = K c , but Wh > and small. According to the 
renormalization scenario discussed above, at some larger length scale the periodic potential 
disappears, and if it was weak initially we end up with an effective theory (15) at the smaller 
cutoff A' with small D(A') and K(A') « 1/4. The recursion relations (16)-(17) then imply 
that at even larger length scales D will renormalize to zero, so the system is in the superfluid 
phase. The same remains to be true at small Wh as long as K < K c {wh) < 1/3, where K c (wh) 
is the initial value of K for which K(A') = 1/3. Since K increases under renormalization 
(i. e. the 0(t> 2 )-term in Eq. (11) is positive), K c < K c (wh) < 1/3 for a small w^. I 
conclude that: 1) there is an intervening Bose glass phase between the superfluid and the 
Mott insulator at weak disorder, and 2) the BG-SF phase boundary must curve upward not 
to intersect K = 1/3 line (see Fig. 1). 

What happens at large disorder Wh can not be concluded on the basis of the perturbative 
renormalization group, since this would eventually correspond to large D in the eqs. (16)- 
(17), which is outside of their domain of validity. I still expect that at any K at arbitrarily 
large disorder the system must be in the localized phase, so the BG-SF phase boundary 
should turn left, as suggested by the dashed line at Fig.l. What happens at smaller K is 
an interesting question, since one is approaching a singular limit of non-interacting particles 
(K = 0). Since in this limit the system should be in localized phase for any disorder, I expect 
the phase boundary eventually to turn downwards, so that there may be a reentrance of the 
localized phase as K decreases at fixed disorder. At the present time however this remains 
a speculation. 

The reader could object that the conclusion that at any disorder the transition from the 
superfluid state is always into a Bose glass may be an artifact of the assumption that the 
disorder distribution is Gaussian and hence unbounded. One can easily show however that 
averaging over a bounded distribution (square shaped, for instance) again gives the -u^-term 
in (9) as the most relevant one, together with terms of higher order in 9. Thus I still expect 
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the dual theory at large length scales to have the same form as in Eq. (15), which leads to 
BG-SF transition. What will change is that a bounded distribution should lead to a finite 
region with the Mott insulator phase, as drawn in Fig. 1, instead of only the half-line at 
Wh = and K > K c , as in the Gaussian case. The transition between the Mott insulator 
and the Bose glass should be first order, since it corresponds to local collapse of the gap [^0| . 

Finally, there is a question of what happens if at the microscopic scale v >> w^. 
At K = K c < 1/4, without disorder one is right on the MI-SF separatrix in the K — v 
renormalization group flow diagram. At small Wh and K = K c it seems plausible than small 
disorder will make v renormalize downward faster than it would without it, so that the 
system would scale towards small v and some K < 1/4. If this is the case, all of the previous 
conclusions would remain the same. Thus even in this case it is natural to expect that the 
transition is again of BG-SF type. 



4 Two dimensions 



4.1 Dual theory 

In 2+1- dimensional lattice theory in Eq. 6 the constraint on M may be resolved by intro- 
ducing another set of integer three component vectors by M = V x N, where again the 
lattice version of the 3D curl operator is assumed. Then: 

Z = lim rn^exp{z27r^/; a -i liQ + |^/J a -5[VxI]}, (18) 

where, 

s = yE( Vx ^- e E^ Vx A + ^ E(v x A)% (19) 

i, a i,a i,a 

and I introduced an auxiliary set of integer three-component vectors {h, a }, the summation 
over which in the limit y — > forces the real valued gauge-fields {^4i )Q J to take strictly integer 
values. To assure the gauge-invariance of the above expression the constraint V • / = must 
be enforced. This is the crucial difference from the ID case, where one could have also 
introduced the two-component vectors {h t0l }, but due to the lack of gauge freedom in 1 + 1 
dimensions, this time without the constraint. So in ID one may sum over {k a } freely and 
end up with the theory (7). In 2+1D however we need new degrees of freedom to enforce 
the constraint. In fixed gauge V • A = the partition function can finally be written as 

Z = lim / Y[[d9 n dA n ] exp{- ]T cos(6 n+fl -0 n - 2nA n ^) - S[V x A]}, (20) 

as can be checked by going backwards from Eq. (20) in Villain approximation and by 
performing the integrations over angles {#„}, which enforces the constraint V • I — 0. Index 
n now labels the sites of a three dimensional lattice, p, — x,y, f, and — n < 9 n < ir. With e 
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finite, one recognizes the last form of the partition function at T = as the 3D "frozen" |p8 
anisotropic lattice superconductor in random magnetic field along one (time) axis. 

The expression (20) is the sought dual representation of the original problem in 2D. 
The angles {6^} represent disorder variables |33|], since they randomize the configurations 



of the gauge-field A. In fact, without the coupling to the disorder variables, the action in 
Eq. (20) would be quadratic in the gauge-field, which may be identified with the sound 
excitation mode in the parent superfluid state of the original bosons. The long-range order 
for ^-variables destroys superfluidity by gapping this sound mode in the pure case. In this 
sense the system described by Eq. (20) is dual to the original Hamiltonian (1): the superfluid 
phase of bosons corresponds to the disordered phase of variables, while the insulating phase 
(MI or BG) of bosons corresponds to the ordered phase of ^-variables. Assuming that by 
leaving small but finite y in the Eq. (20) the universality class of the transition at T = is 



not changed [2J|, 24], one may write the soft-spin version of the lattice gauge theory in (20) 
Ml: 



as 



S= Jd 2 fdz[\(V-tqA(f,z))^(f,z)\ 2 + /j 2 \^(f,z)\ 2 + -\^(f,z)\ i + (21) 

^(V x A(f,z)) 2 + |(V x A(r,z))j + h(r)(V x A(r,z)) s ], 

which is a Ginzburg-Landau theory for a strongly type-II superconductor (assuming A >> q 2 ) 
in transverse gauge V ■ A = 0, in the random magnetic field along the z-axis. Non-zero 
parameter 7 indicates the inherent anisotropy of the quantum problem. The transition is 
tuned by varying the mass /x, and the condensation of the complex disorder field \I> signals 
the destruction of superfluidity. 



4.2 Renormalization 

In the field theory (21) disorder appears as the random magnetic field along the ^-direction. 
This term now plays the role similar to the ttvterm in ID. Note that the random magnetic 
field depends only on two (spatial) coordinates and is completely correlated in time, as typical 
for quantum problems with static disorder P5|] . The quantity of interest is the disorder- 
averaged zero-temperature free energy F = —InZ. Using the replica trick to average over 
disorder, one arrives at the action for interacting replicas of the system 

S' = E / d 2 rdz{\{V - iqA a )m a \ 2 + /i 2 |vig 2 + ^|* Q | 4 + -(V x A a ) 2 + 

a=l ~~ " 




Similarly to ID, the effect of averaging over disorder is to introduce an off-diagonal (in replica 
indices) contribution to the gauge-field propagator, as determined by the last term in Eq. 
(22). 
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The point of transition in the dual theory is reached by tuning the renormalized mass 
of the disorder field m to zero. Simple power counting tells that one needs to consider four 
dimensionless coupling constants: the dual charge q 2 = q 2 /p, the quartic term coupling A = 
X/p, the anisotropy parameter 7 and the width of the random field distribution Wh = Wh/p, 
where p is an arbitrary infrared scale. However, little analysis shows that this is not enough, 
and to complete the theory one needs to include one more coupling constant. Consider the 
diagram in Fig. 2: since the gauge-field propagator has off-diagonal part in replica indices, 
this diagram generates the quartic term which couples different replicas even though this 
term is absent initially. This is analogous to the way the standard quartic term A would get 
generated by the fluctuations of the gauge-field in the pure theory. We already saw that a 
similar generation of a new disorder term happens in ID theory as well. To complete the 
theory one therefore must add the following term to the action in Eq. (22): 

S " = 4t/ d 2 rdzdz'\^ a (f,z)\ 2 \^^r,z')\ 2 , (23) 

a,P=l J 

with w > 0. Note that this term has precisely the form which would arise from averaging 
over an independently random mass in (21), correlated along the z-direction and with the 
Gaussian distribution with a width w. Thus the full theory S = S' + S" contains five coupling 
constant at the critical surface m = 0: the fifth dimensionless coupling is w = w /up 2 . 

To address the quantum (T = 0) critical behavior at the SI transition I consider the 
renormalized theory right at the transition point m = 0: 

N , w / N 

1 T i/i x — ■» / /— \iOi t /— /\|2 



Sr = E / d 2 rdz{Z*\Z^V» - iqA a> ^ a \ 2 + -\^ a \ A - ^- E / dz'\V a (r, z)\ 2 \* p (r, z 1 

+^(V x A a ) 2 + ^(V X JL)]-^EJ dz,{ y x Mr,z)),(V x Ap(f,z%}m 



The lowest order contributions to the renormalized quantities are determined by the diagrams 
at Fig. 3. The infrared divergences are regulated by evaluating quartic vertices at the usual 
symmetric point 

ki ■ kj = '—^ , (25) 

with an additional condition k\ )Z = k<i )Z = —k^ z = —k^ z , necessary because the disorder is 
correlated along the z-axis, so that individual k z components, and not only their sum, are 
conserved in w- vertices. The polarization is evaluated at a finite momentum c • p, where c is 
a constant to be specified shortly. The result of the one-loop renormalizations is: 

Z* = l- hf (26) 
Z z = 1 + \w (27) 
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w' 



Za = 1 + ^q 2 (28) 

A' = A — h^tlx 2 - -^g 4 + 2Xw(V2 ln(l + y/2) + 4= ln (V2 + v^)), (29) 
o 2y2 v6 

«; + 2^ 2 (-L ln(l + y/2) + ln(>/2 + v^)) - ^At&, (30) 

Z~i = Zh = Z x = Zy = 1. (31) 

The non-trivial numerical factors in Eqs. (29)- (30) arise from integration performed directly 
in 3D and from the choice of the renormalization point (25). 

Let me first discuss the renormalization of the anisotropy parameter 7 and the random- 
field disorder variable w^. The renormalized couplings 7 r = ^{Z 1 /Za and xbh >r = xbhZ h /ZA 
obey the equations: 

^ = -7,(^ 2 + 0(g 4 )) + 0( 7r 2 ), (32) 
d ^ = w h Al-VA) + 0(wl), (33) 

where t = — ln(p) and 

VA = (34) 

is the anomalous dimension of the gauge-field propagator. A small anisotropy appears to be 
irrelevant at any fixed point with a finite dual charge. I checked that this remains true when 
the term 0(q 4 ) is included in Eq. (32). Note that anisotropy in the gauge-theory can not 
simply be rescaled out, as it would be possible in the XY model. Although the dual theory 
(21) should be strongly and not weakly anisotropic in the continuum limit a — > 0, e —>■ 0, 
I will assume that there are no stable anisotropic (7 7^ 0) fixed points in the theory. This 
would be in agreement with what is known for the classical scalar electrodynamics close to 



four dimensions for large number of complex field components |36| . I therefore set anisotropy 
7 = hereafter. 

The scaling equation for the renormalized dual charge g 2 = q 2 jZ& may be written 
exactly as 

= i 2 (l-T7A), (35) 

as follows from the definition of the anomalous dimension of the gauge- field propagator. 
Thus at any fixed point in the theory with a finite dual charge, tja = 1 [ 57|| - In particular, 



this implies exact vanishing of the linear term in Eq. (33) at the MI-SF fixed point of the 
pure system. The fact that the scaling dimension of the random-field disorder Wh is simply 
1 — f]A follows from the observation that the transverse part of polarization is diagonal in 
replica indices to the lowest order in w^. The term in Fig. 4, for example, does not contribute 
to Z h since it turns out to be completely longitudinal. 

The fate of the random-field disorder vbh is determined by higher order terms in the 
Eq. (33). Marginality of Wh follows from the quantum nature of disorder: if disorder was 
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correlated along slightly less than one dimension its scaling dimension would have been 
negative. Guided by the ID problem, I will make a simplifying assumption that apart from 
generating the quartic term w, disorder coupling Wh is irrelevant at the MI-SF fixed point. If 
this assumption turns out to be incorrect, the calculation presented here would be pertinent 
to the crossover regime towards a new critical point with Wh 7^ 0. However, since Wh can at 
worst be marginally relevant at the pure fixed point, and therefore can initially grow much 
slower than w, we are assured of a region for a small initial tih where ignoring its possible 
growth is a sensible approximation. 

The scenario in 2D is thus similar to what we already encountered in ID: initially 
disorder in the problem is represented by the random magnetic field Wh, which I assume is 
marginally irrelevant at the MI-SF fixed point. Before renormalizing to zero Wh generates 
another disorder-like coupling w which is relevant and which grows. In the rest of the article 
I will then set Wh = 0, and follow the renormalization of the remaining three coupling 
constants. 

The one loop /9-functions for the remaining coupling constants g 2 , w r = w' jZ\ and 
A r = W jZ\ are: 

-dt= q r~ 16^' (36) 

^ = A r (l + \ql + 2(V21n(l + V2) + ± ln(V2 + V3))w r ) - ^|±\ 2 - ^=g r 4 , (37) 

^ = w r (2 + l -ql - ^j=\ r ) + 2(i= ln(l + y/2) + -±= m(V2 + y/3))w 2 r . (38) 

For w r = these equations reduce to those studied previously in the context of critical 
behavior of superconductors [[37]]. In that case for a choice c > 5.17 there are four fixed points 



in the theory: Gaussian (g 2 = A r = 0) and 3D XY (g 2 = 0, \ r ^ xy = 2.09), both unstable 
in direction of the dual charge, tricritical (g 2 = 16/c, A r = A_), unstable in A-direction, 
and the MI-SF critical point (g 2 = 16/c, A r = A + ), believed to be of "inverted" XY type 
p9| (see the inset on Fig. 5). The Gaussian and the tricritical fixed points are connected 
by the straight separatrix which determines the tricritical initial value of the Ginzburg- 



Landau parameter k = yA/2g 2 = k c in the superconducting problem. The tricritical value 
of k has been estimated both analytically |38| and numerically f39|, and we may tune the 
renormalization point to match that number: k c = 0.42/a/2 requires c = 20. It is worth 
noting that for w = our one-loop approximation gives encouragingly good estimates of the 
critical exponents both at the unstable XY [y = 0.63) and at the stable, charged fixed point 
(z/ = 0.61, r] = -0.20) [0. 

A small disorder w is relevant at the MI-SF fixed point, as follows from the last (3- 
function, and more generally may be expected on the basis of Harris criterion |^TJ. Both 



perturbative |37| and non-perturbative pq| , |29|] , |30l calculations of the correlation length 
exponent at the pure critical point yield values v < 1, suggesting relevance of the correlated 
random-mass disorder. Starting from the strongly type-II region (A >> g 2 ), the flows are 
attracted by the random critical point at g 2 = 0.8, w c = 3.71 and A c = 29.72, which I identify 
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as BG-SF critical point (Figure 5). The BG-SF critical point exists for any choice of the 
parameterization by c, unlike the charged fixed points in the pure theory. This is similar to 



the disordered classical scalar electrodynamics close to 4D [42], where randomness restores 
the critical point in the theory. The flow towards BG-SF critical point is oscillatory, which 
should lead to oscillatory corrections to scaling. 

4.3 Critical exponents 

Having obtained the approximate renormalization group flow and the fixed points, I proceed 
to calculate the critical exponents. First, note that the correlation length exponent v and the 
dynamical exponent z may be assigned a purely thermodynamical meaning, since the BG-SF 
transition temperature T c ~ 5 ZU and the zero-temperature compressibility k ~ § v \ d ~ z ) . 5 is a 
parameter measuring the distance from the BG-SF transition at T = 0. One may therefore 
calculate these exponents from any representation of the partition function. The correlation 
length critical exponent v is naturally defined away from the critical surface m — 0, while 
the calculation here is performed precisely at the critical surface. Nevertheless, I may still 
obtain v from the knowledge of yet another renormalization factor Z^i [25] that accounts 



for the renormalization of the |\1/| 2 term in the theory (22), evaluated at the critical surface 
m = 0. I find: 

= 1 + \ ~ f , (39) 

to the lowest order in the coupling constants (Figure 3 a). In analogy with the thermal 
critical phenomena, one may define the exponent 7 as: 

7 = l-lim— — §-, 40 

so that the exponent v then follows from the standard scaling relation: 

v = (41) 
2-Vd 1 ^ 

rjd is the anomalous dimension of the vortex field \1/ 

% = -lim— -A 42 

p^o amp 

Using the lowest order eqs. (26) and (39), expressing the bare coupling constants in terms of 
the renormalized ones, and taking the infrared limit p — >• 0, the correlation length exponent 
becomes a series in fixed point values of the renormalized coupling constants. To the lowest 
order I obtain: 

,= i(l + ^4^) = 1.38. (43) 
Similarly, the dynamical critical exponent z is defined by: 

z* 1 = 1 + lim — — ^. 44 
p-o dlnp 
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From Eq. (27), to the lowest order in the critical coupling constants I find: 

z = l + -Z = 1.93. (45) 
4 

At the BG-SF critical point v > 1, as expected ^T|. The dynamical exponent is very 
close to two, which was conjectured to be exact in 2D, based on the assumption that 
compressibility at the BG-SF transition stays finite. The values are also in good agreement 
with the Monte Carlo results of ref. JT?j [y = 0.9 ± 0.1 and z = 2 ± 0.1) and the real-space 
study of ref. @ (v = 1.4 and z = 1.7). 

In contradistinction to the critical exponents v and z, the standard exponent rj char- 
acterizes the behavior of the correlation function for the superfluid order parameter at the 
transition, and can not be straightforwardly obtained from the dual theory. 



4.4 Universal conductivity 

Since the calculation is performed right at the transition point (where the mass of the dual 
field m = 0), one may calculate the universal dc conductivity in similar way as the critical 
exponents. Before proceeding however, two important points need to be clarified. First, 
the dual theory describes the dirty-boson system at T = 0, and therefore the universal 
conductivity I will calculate corresponds to T — > 0, ui — > order of limits. The reader should 
note that the order of limits in typical experiments is reversed, which raises subtle issues 
about the actual dissipation mechanism tested by the experiments to which I will return 
in the concluding section. Nevertheless, one may compute T = critical conductivity as a 
fundamental characteristic of the BG-SF fixed point. The second point concerns my intention 
to obtain the conductivity of bosons from the dual theory for the vortex field. This question 
is intimately linked with the mechanism of dissipation in the quantum (T = 0) problem. 
As already mentioned, the dual theory describes the disappearance of superfluidity due to 
destruction of the gapless mode, which, as argued, can happen only through condensation 
of vortices. The physical picture is that on the insulating side of the transition vortex loops 
in 3D classical electrodynamics blow up, completely disordering the phase of the superfluid 
order parameter. The ground state of the Bose glass insulator is therefore a superfluid of free 
vortices and antivortices, which represent the 2D space projection of the 2+1D vortex loops, 
while the original bosons are localized by random potential. In contrast, in the superfluid 
ground state bosons are phase coherent while vortices and antivortices are tightly bound 
into pairs. Right at the transition it seems natural to assume that both bosons and vortices 
should be mobile which then leads to a finite metallic conductivity [|Kj. Consider the current 
of bosons 

h = e*N b (46) 

from one end of the sample to another; motion of vortices perpendicular to the direction of 
the boson current induces the voltage by Josephson relation: 

V^ k A (47) 
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where N v is the vortex flux. The boson resistance is therefore 



Rb 



N v h 



(4£ 



On the other hand, vortices "see" bosons exactly the same way as bosons see them |[44|| ; 
writing the above relations for current and "voltage" of vortices instead of bosons we find 
that: 

V v N b h 



R v 



(49) 



where e v is an arbitrary " charge" of a vortex. Combining the last two equations I find 



h 



el/ h el 



(50) 



where cx fo and a v are boson and vortex conductivities, respectively, and I used the fact that 
conductivity equals conductance in two dimensions. 

The last expression enables one to calculate resistivity of the bosons a^ 1 by using the 
Kubo formula for vortex conductivity a v . The derivation of the Kubo formula in the theory 
with replicas is presented in Appendix B. The simplest way to evaluate the integrals over 
the correlation functions in Eq. (71) is to use the scheme of dimensional regularization. 
Recalling that the couplings A and q 2 have dimensions of frequency, it follows that in the 
static limit ou — > to the lowest order in these coupling constants we need to calculate the 
integrals in Eq. (71) only at k = 0, so that they will give no contribution to a v to this order. 
To the lowest order a v can then be written as (k = (0, oo)): 



<j„ 



a v n + a 



V,lj 



(51) 



where 



lim — ( 



d?q 1 
(2tt)V 



(27T) 2 q 2^ + ky' 16 



8 f d 3 q S(g) 



Aw f d 3 qd 2 p± 



q x Px 



in dimensional regularization, in agreement with the result of Cha et al. 1 10 
by the diagrams at Fig. 6: 



(2tt) 5 q*(jfi_ + <£)($+ k) 2 (j? ± + (q, + u) 



(52) 

. a %1 is given 
(53) 



in the static limit oo — > 0. The second term in a v> i vanishes due to q x integration. The 
self-energy appearing in the first term is to be calculated to the lowest order in w at the 
transition point m = 0. For small momenta it can be written as 



w 
2~ix 



In | & 



(54) 
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as follows from the result for the dynamical exponent z. Inserting the last expression into 
Eq. (53) and introducing Feynman parameters to integrate over q x and q y leads to 



^ = 4^" (55) 



The integral in the numerator is: 



I = [°° d q z f 1 d A 1)ln|gz| = 5.89. (56) 
i-oo H Jo q 2 z + 2tq z + t V ; 

Recognizing the combination w/nu 2 as the dimensionless coupling w, the lowest order result 
for the dc vortex conductivity becomes: 

ff " = <i + ^r4 (57 > 

and therefore formally divergent due to the infrared (static) limit. Rewriting a v as a series 
in the renormalized coupling w r however, to the lowest order 

* v = Z^£ + ^\imw r ) e j, (58) 

where I introduced the wave-function renormalization factor Zq, to account for the fact that 
w r describes correlation functions of the rescaled vortex field. To the lowest order the boson 
dc resistivity is then determined by the critical coupling constants as: 

err 1 = (-H q 2 A w c ) — . (59) 

Numerical values of the critical couplings as found from the one-loop /9-functions eqs. (36)- 
(38) yield the estimate of the universal conductivity of bosons at the BG-SF transition in 
2D: 

e 2 

a b = 0.25-i. (60) 



5 Summary and discussion 

In the present paper it was shown that the system of disordered interacting bosons at a 
commensurate density has a superfluid-insulator transition at zero temperature governed by 
the random fixed point, which was interpreted as a sign of a BG-SF transition. This critical 
point can be naturally identified within the dual theory for the transition which describes 
topological defects in the ground state. The coupling constant which represents the disorder 
relevant at the pure MI-SF fixed point is not initially present in the microscopic dual theory, 
but becomes generated at larger length scales by the microscopic disorder and the lattice 
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potential. I discussed the phase diagram in ID, and calculated approximately the critical 
exponents v and z and the universal conductivity at the BG-SF transition in 2D. 

I argued that only in the dual theory the BG-SF critical point may become apparent. 
It is not obvious to which terms would the generated disorder in Eq. (13) correspond in 
the theory written in terms of the superfluid order parameter, for example. The situation 
is reminiscent of the Kosterlitz-Thouless transition in the 2D XY model: while in principle 
the information on the transition is contained in the Ginzburg-Landau- Wilson functional, 
its extraction is hopelessly difficult; on the other hand, by focusing directly on the behavior 
of vortices the problem of the critical behavior can be solved exactly. It appears the same 
is accomplished by the dual formulation of the dirty boson problem, at least in ID where 
the random critical point is perturbatively accessible. In 2D the dual theory has a strong- 
coupling critical point, and one can obtain the critical quantities only approximately. 

It is interesting how the relevant disorder coupling is generated by the random chemical 
potential and the lattice, although they are both ultimately irrelevant at the BG-SF fixed 
point. In ID the recursion relations imply that weak lattice potential always renormalizes 
to zero, but that while doing so it conspires with random chemical potential to induce 
another disorder term in the theory. At some larger length-scale the initial disorder becomes 
redundant in ID, and one is left with an effective theory in the familiar Haldane's form. In 
a similar way, in 2D the random magnetic field generates the random mass for the vortex 
field, which immediately becomes the most relevant perturbation at the pure fixed point and 
whose flow ultimately determines the critical behavior. I assumed, but could not prove, that 
in 2D the initial disorder in form of the random magnetic field also becomes irrelevant once 
it generated the important random mass. 

The proposed phase diagram in ID has one counterintuitive feature: there is a region for 
K > K c where by increasing disorder one turns a Bose glass insulator into a superfluid. This 
comes as a consequence of renormalization in the theory and I believe is a genuine effect. The 
phase diagram in Fig. 1 is in disagreement with the density-matrix renormalization group 
study of Pai et al. []13"| , who claimed the direct MI-SF transition at small disorder. This 
work has been criticized recently by Prokof'ev and Svistunov |45] who performed a Monte 
Carlo calculation to find the phase diagram in qualitative agreement with the one discussed 
here. 

The performed renormalization group calculation of the critical exponents in 2D, al- 
though still approximate, has a slight advantage when compared with other ad hoc schemes 
in that it leaves room for a systematic improvement by going to higher orders in perturbation 
theory. Although this is possible in principle, it is rather complicated since there are three 
couplings to consider. The relatively good agreement between my results and those of other 
approaches is therefore somewhat fortuitous. The reader should also note that the values of 
the critical exponents are close to the results of other studies [17|, [19| which are performed 
at incommensurate densities of bosons. In ID it is easy to show that nothing in the conclu- 
sions is changed by allowing h ^ 0, since by shifting the dual field one may always place h 
under the cosine term in Eq. (9), so it disappears together with the lattice potential at large 
length scales. In 2D incommensurability corresponds to non-zero magnetic field along one 
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direction, the effect of which at the random fixed point is presently not clear. Physically one 
would expect the incommensurability should become irrelevant at the BG-SF critical point, 



although it would be desirable to have a concrete calculation in support of this claim [46 . 

The universal conductivity at the BG-SF critical point in Eq. (60) turns out to be 
somewhat larger than in the Monte Carlo calculations of Wallin et al. ]T7|] , and smaller than 



in Batrouni et al |16[ . It is smaller that at the MI-SF transition |10| , ||11|| , which seems 
plausible. It is also smaller than what is found in most experiments, where a ~ (2e) 2 /h, and 
it has recently been proposed that this discrepancy may not be a coincidence. Damle and 
Sachdev |^3| pointed out that at the critical point a becomes a non-trivial universal function 
of ou/T, and that the experimentally relevant limit corresponds to behavior of this function 
at small values of its argument. In this regime the relevant dissipation mechanism is the 
random scattering of thermally created quasiparticles, rather that the coherent transport of 
the field-induced defects in the ground state, which would correspond to the opposite limit 
uj/T — > oo. The direct comparison of the conductivity calculated here with the present 
experimental data is therefore probably not appropriate. I may still note that number in 
Eq. (60) is intriguingly close to recent measurement in thin Pb film f471| : o c « 0.25(2e) 2 //i. 
The experiment which would measure the high-frequency conductivity at low temperatures 
would be closer to probing the quantum regime for which the result in Eq. (60) is directly 
relevant. 

Finally, the kinship between the BG-SF critical points in ID and 2D suggests that one 
could try to approach the strong-coupling critical point in 2D via an expansion around ID 
where the BG-SF fixed point happens to be located at zero disorder. Although the apparent 
dissimilarity between the forms of dual theories (9) in ID and (22) in 2D at first sight 
makes this idea unlikely to work, all that is needed is analyticity of the critical exponents 
as functions of the dimension for 1 < D < 2, for such an expansion not to be impossible in 



principle. Some steps in this direction have already been taken in [4"5[], and this will be a 
subject of a future publication (49| . 

This work has been supported by NSERC of Canada and Izaak Walton Killam founda- 
tion. 



6 Appendix A: dual field propagator in ID 

From the dual theory (9) the inverse of the propagator for the field 9 is: 

&- 1 = K(u 2 + c 2 k 2 )I -wk 2 5{uj)M, (61) 

where G is N x N matrix, / is a unit matrix, and M is also a. N x N matrix with all elements 
equal to one. The second term is due to disorder, and the delta-function in frequency 
originates in the quantum nature of the problem: disorder is represented by a potential 
random in space, but completely correlated in imaginary time. By noticing that M 2 = NM 
we may write 

^ift^W^S^"^- (62) 
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In the relevant limit of vanishing number of replicas N — > the propagator therefore takes 
a simple form: 

G = Ki^TdWj 1 + lw M ' (63) 

7 Appendix B: Kubo formula for disordered dual the- 
ory 

We want to find the response of the system described by Eq. (22) to the external vector 
potential a which couples to the disorder field via minimal coupling. I work in units where 
the charge of disorder field e v which determines the strength of the coupling to the external 
probe (not to be confused with the running dual charge q) and H are both set to unity. With 
the external probe, the kinetic part of the action in Eq. (22) becomes: 

N 

S ' = E / d 2 fdz\(\/ - iqA a + ia)^ a \ 2 , (64) 

with the rest of the action unchanged. Disorder-averaged current of the disorder field is then 
given by: 



(j(r,z)) = - \ s=0 . (65) 

oa{r, z) 

Angular brackets denote quantum average over the ground state and overline stands for the 
average over disorder configurations. The vortex conductivity is defined by linear response 
to the vortex "electric field" as: 

(Mz)} = -J d 2 7dz'a v (r- ?, z - z') dn f d f\ (66) 

Integrating by parts the last equation and going to Fourier space I find: 

/ \ 1 f ,2^, /• x 5 2 F\a\ . „. 

a v (uj) = d rdzex.p(iuz) \ s=0 , (67) 

uj J 5a x (r,z)da x (0,0) 

where I chose a to be along x direction. Using the replica trick to write F and performing 
the differentiations in (67) I obtain: 

^(u;) = - lim 4{2£ / exp(- 1 S)|^ a (0,0)| 2 -/ / d 2 fdzexp(tujz) exp(-S) 3x (f, z)j x {6, 0)} 

(68) 

where the current is given by the standard expression summed over replicas: 

N 

3 = E(**aVtf Q - + 2q\^ a \ 2 A a ). (69) 
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Since 

= hfE^A = [ eM S)A, (70) 

in the limit N — > 0, inserting the above expression for the current, the last equation can be 
rewritten in the final form: 

~^ V2 E d%d3 f f tlf d¥2 ^a(gl)^(g2)^(pi)^(P2)A a , a (fe + & - ? 2 )A^(pl -ft - £))},(71) 

with /c = (0, a;), V is the 3D volume, and I restored the correct units. Apart from appearance 
of the replica indices, first two terms are the same as in Cha et al. |]10| . Simple power-counting 
in the last expression implies that o v does not scale with frequency, and is thus given by 
an universal constant, in its natural units e^/h. This conclusion remains true in the fully 
renormalized theory as well, since the current operator can not acquire anomalous dimension 
due to gauge invariance. 
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Figure Captions: 



Figure 1: The phase diagram for the commensurate ID Bose-Hubbard model in random 
chemical potential. K = JU/t/2'K 2 and Wh measures the strength of disorder. Dashed line 



represents a more speculative part of the BG-SF phase boundary. The BG-MI transition is 
expected to be first order. For further discussion see the text. 



Figure 2: Diagram of the order w^q 4 which generates the quartic interaction w between 
different replicas in 2D. Dashed line represents the gauge field propagator. 



Figure 3: The lowest order contributions to the self energy a), A' b), w' c) and Za d). Dashed 
line represents the gauge field propagator, the wavy line stands for the interaction vertex w 
and dots denote the interaction vertex A. 



Figure 4: Purely longitudinal contribution to the off-diagonal (in replica indices) polarization. 



Figure 5: Renormalization group flow and the fixed points in q 2 = q 2 plane. The inset shows 



the flow in w = plane |j3~T |. 



Figure 6: Contributions of the order w to the vortex conductivity a v 
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